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I. INTRODUCTION 



In the limit of infinite heavy quark masses, the weak decay form factors describ- 
ing semileptonic transitions between any two ground-state (pseudoscalar or vector) heavy 
mesons, M[v) —>■ M'{v')ii', are described by a universal form factor ^(y). This so-called 
Isgur-Wise function depends on the velocity transfer y = v ■ v' and is normalized at zero re- 
coil, ^(1) = 1, where the initial and final meson have the same velocity [Q. These remarkable 
results follow from an implicit spin-flavor symmetry, which QCD reveals for heavy quarks 
although it is not explicit from its Lagrangian . An expansion about this symmetry limit 
is afforded by the construction of the heavy quark effective theory (HQET) In the 

infinite mass limit, its effective Lagrangian is explicitly invariant under spin- flavor symmetry 
transformations. HQET thus provides a convenient framework in which to analyze the prop- 
erties of hadrons containing a heavy quark. In particular, it allows a systematic expansion 
of weak decay form factors in powers of 1/ uiq • 

To leading order in this expansion one recovers the Isgur-Wise limit, in which a large 
set of otherwise unrelated form factors reduces to the Isgur-Wise function. This function 
describes the dynamical properties of the cloud of light quarks and gluons surrounding the 
static heavy quarks. Being a hadronic form factor, it can only be investigated using non- 
perturbative methods. One such method is provided by QCD sum rules 0, which were 
originally developed for light quark systems and have yielded many nice results which are 
competitive with lattice computations. Recently, several authors have used QCD sum rules 
to calculate hadronic matrix elements in HQET P, p!^JTT| , p!2|Jl3| , p^JT5| , [l6|JT7| , p!8[| . The more 
refined of these analyses included radiative corrections to both the perturbative and non- 
perturbative contributions. For the Isgur-Wise function, however, the radiative corrections 
to the perturbative part of the sum rule were only incorporated in leading logarithmic ap- 
proximation the complete two-loop corrections to the triangle quark loop were never 
calculated for the case of heavy quarks. On the other hand, from the well-studied case of 
pseudoscalar decay constants it is known that next-to-leading logarithmic corrections can 



be quite substantial and should in principle be taken into account [11,12,13 



With the advance of HQET, considerable progress has been made in the calculation of 
radiative corrections. In Ref. |l^, for the first time an exact two- loop result was obtained 
for a heavy meson form factor, in this case for one of the universal functions that appear 



at order 1 / tuq in the heavy quark expansion. The authors of Ref. developed a general 
method to compute the first two terms in an expansion of a two-loop diagram in HQET as 
a power series in (y — 1), and applied their technique to obtain an expansion of the two-loop 
corrections to the perturbative part of the Isgur-Wise function close to zero recoil. 

In this paper we derive the exact result for the two-loop corrections to C,{y)- To this end, 
we develop several techniques to evaluate two-loop integrals in HQET involving two heavy 
quarks with different velocities {v and v') and different residual momenta {k and k'), among 
them the method of integration by parts ||T9| , |20| , pl| and the use of differential equations 



We also introduce an integral representation for a general two-loop diagram which is 
particularly convenient for QCD sum rule calculations. These techniques are rather general 
and can be readily applied to other cases. In Sec. |1| we briefly review the sum rule analysis 
of the Isgur Wise function. The calculation of the two- loop perturbative corrections to ^ (y) 
are described in Sec. fT|. We discuss in detail the contribution of each individual diagram, 
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so that the interested reader can follow the analysis step by step. After renormalization, we 

m 



compare our exact result with the expansion around zero recoil given in Ref. 



and find 



agreement. Sec. |V| deals with the renormalization- group improvement and the numerical 
analysis of the sum rule. We find that the effects of radiative corrections to the Isgur-Wise 
functions are moderate and well under control. Sec. contains the conclusions. 



II. SUM RULE FOR THE ISGUR-WISE FUNCTION 

The derivation of the QCD sum rule for the Isgur-Wise function has been dealt with 



at length in Refs. [|10],|12|,|14| , to which we refer the interested reader for details. Here we 
restrict ourselves to a brief review for the purpose of introducing the necessary notations 
and recalling the main ideas of the method. One studies the analytic properties of the 
three-current correlator 

J dx dz e^e^'-^-) {0\T{[qTM'h'l,[h'Th]^,[hTMql}\0) 

= E{uj,uj',vv')Tt{V'TV}, (1) 

where h and h' describe heavy quarks with velocities v and v' in the effective theory. These 
quarks have "residual" momenta k and k', which are related to the total external momenta 
by P = rnqv + k and P' = rriQ/v' + k', where mg and rriQ/ are the heavy quark masses. 
Depending on the choice Tm = —75 or = 7^ — v^, the heavy-light currents interpolate 
pseudoscalar or vector mesons, respectively. The Dirac structure F of the heavy-heavy 
current is arbitrary. Usually, however, this is a flavor-changing weak current, in which case 
F = 7^(1 — 75). The Dirac structure of the correlator is entirely contained in the trace over 
"spin wave functions" 

P.i±^r„. r^r„,i±^. (2) 

which act as projection operators: 

fV = V = -Vf, 
-^'p' = p' = py. (3) 

The coefficient function S(a;, uj' ,v-v') in (|l]) is analytic in the "off-shell energies" uj = 2v-k 
and uj' = 2v' ■ k', with discontinuities for positive values of these variables. In particular, 
it receives a double-pole contribution from the ground-state mesons M and M' associated 
with the heavy-light currents. This pole is located at u = u' = 2A, where A = rriM — rnq = 



rriM' — raqi. The residue is proportional to the Isgur-Wise function. It follows that [|12 



[ijj — 2K + ie)[ijj' — 2A -I- le) 

where y = v ■ v', and F corresponds to the scaled meson decay constant in the effective 
theory {F ~ fMy/f^M)- Note that both F and the Isgur-Wise function are defined in terms 
of matrix elements in the effective theory and are therefore scale-dependent quantities. 
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In the deep Euclidean region, the correlator can be calculated perturbatively by using the 
Feynman rules of the heavy quark effective theory The idea of QCD sum rules is that, 
at the transition from the perturbative to the nonperturbative regime, confinement effects 
can be accounted for by including the leading power corrections in the operator product 
expansion of the three-point function. They are proportional to vacuum expectation values 
of local quark-gluon operators, the so-called condensates [||. One then writes the theoretical 
expression for the correlator in terms of a double dispersion integral, 

Sth(cf , uj\ y) = [du du' ^^'^^ ' — — + subtractions, (5) 

J [u — uj — ie)[u' — to' — ie) 

and performs a Borel transformation in u and u' (see Appendix A for the definition of the 
Borel operator). This yields an exponential damping factor in the dispersion integral and 
eliminates possible subtraction terms. Because of the flavor symmetry it is natural to set 
the associated Borel parameters equal: r = r' = 2T. Following Refs. |]I2|,|T5[, one then 
introduces new variables u!± = |(i/±z/'), performs the integral over u^, and employs quark- 
hadron duality to equate the integral over up to a threshold cuq to the Borel transform 
of the pole contribution in (^). This yields the Borel sum rule 

F\fi) e-^^/^ = J d^+ e-'^+/^ p,h(a;+, y) = K{T, uo, y) . (6) 



The effective spectral density pth arises after integration of the double spectral density over 

OJ-. 

To lowest order in perturbation theory, the theoretical expression for the right-hand side 
of the sum rule is given by 



/ 9 \ 2 7 

K{T^ ^0, = ^ [--^ j du^ ul e--l^ - (m) 





/ y-l \ {asGG) (2|/ + 1) mg (gg) 

\y+lj ASttT 3 4T2 ■ ^ ^ 

We have included the leading nonperturbative contributions in the operator product ex- 
pansion, which are proportional to the quark condensate (dimension d = 3), the gluon 
condensate {d = 4), and the mixed quark-gluon condensate {d = 5). In the numerical 
analysis in Sec. |^ we will use the standard values (at /i = 1 GeV) 

{qq) = -(0.23GeV)^ 
(asGG) = 0.04 GeV^ 
(gsqa^^G'^'q) = ml {qq) , = 0.8 GeV' . (8) 

At zero recoil, a Ward identity relates the three-point function (|l|) to the correlator of 
two heavy-light currents, from which one derives the sum rule for the parameter F. It allows 
one to replace the product F2g-2A/T ^ /^(^j^ ^jq, 1). Then the final form of the sum 
rule for the Isgur-Wise function explicitly reveals its normalization at zero recoil: 
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(9) 



In the following section we will derive the complete expression for the perturbative cor- 
rections to the function K{T,uJo,y) arising at order as- The one-loop corrections to the 
quark condensate were calculated in Ref . [0 . The mixed and gluon condensate are already 



of order gs or g'^, and consequently one does not have to include radiative corrections to these 
terms at order as- What is missing are thus the order-a^ corrections to the perturbative 
contribution. There are restrictive constraints on the result of this two-loop calculation. 
The normalization of the Isgur-Wise function requires that, at zero recoil, one must recover 
the expression for the two-loop corrections to the perturbative part of the sum rule for F- 
This implies [p!l| , p!^ ,p!^ 



as 

7T 



■ fi 47r2 17 
2 In + — + — 
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(10) 



Furthermore, the two-loop calculation must reproduce the known anomalous dimension of 
both F{fi) and C,{y,fi) • We will see at the end of Sec. |T| how these constraints are 

fulfilled. 



III. TWO-LOOP CALCULATION 



The two-loop corrections to the perturbative contribution to the three-current correlator 
(0) are shown in Fig. |l|. We shall analyze these diagrams separately below. Throughout the 
calculation we use Feynman gauge. For practical purposes, it is useful to realize that the 
dependence of the perturbative spectral density ppert{^+,y) on is known on dimensional 
grounds: 



ppert(t^+>l/) = 



da;+e--+/^Ppert(^+,y) = 2T^ 



pi{y) + p2{y) In — 



pi{y) + p2{y)[ in^ + 7s-- 



(11) 



The coefficient functions Pi{y) are independent of u^. It thus suffices to calculate directly 
the Borel transform of the correlator, corresponding to the second line in this equation. 
The spectral density can then be read off immediately. Below, we will always denote Borel- 
transformed quantities by a "hat". 



A. Gluons attached to heavy quark lines 

Consider the three diagrams Di to in Fig. |I|. The evaluation of the first graph gives 
Di = 16iV,C^(?2^Tr{p'rP7"} 

(12) 



X / ds dt 



{to + 2v- s){uj + 2v ■ t){uj' + 2v' ■ s){uj' + 2v' ■ t) {s - tf 
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where Cp = {N^ — l)/2Nc, and ds = {2'K)~^d^s. The two-loop integral is most conveniently 
performed by using a Fourier representation for the light quark and gluon propagators, 
and an exponential integral representation for the heavy quark propagators. A detailed 
description of this method, as well as its application to the above integral, can be found in 
Appendix A. After Borel transformation we find 



[2(y + l)r/^-^G(0,0,f-l;,) 



(D-A) 



where we have abbreviated 

_ IQN.Cpg'^s (2T)^^-5 T(D)T(D UTrlv'vV] 
^- ^A^)D [2(y + l)r^^^^^^^^ l)Tr|PrP}, 

Here and in the following we will often encounter integrals of the form 



(13) 



(14) 



du ■ 



1 + 2yu + u^) 



(15) 



which are related to generalized hypergeometric functions of the velocity transfer y. 
The contribution of the second diagram is 

D2 = 16N,Cf gl Tr{ V'V V 7°} 
X dsdi 



{uj + 2v s)2(cu + 2v ■ t){uj' + 2v' ■ s) {s - ty ' 



(16) 



Introducing a shifted loop momentum hj t' = t + s, one can first perform the integral over 
t' and then carry out the integral over s by using the master equations ( |A1| ) and (^3) for 
one-loop integrals in HQET (see Appendix A). After Borel transformation, the result is 



Do 



2A 



{D-A){D-3) 



[2(y + 1)] 



D/2-2 



(17) 



Obviously the third diagram, D^, gives the same contribution. 

Next we set D = 4 + 2e and expand in e, using some of the integrals collected in 
Appendix B. We find 



1=1 



1 



J2D, = A\- 2-yr{y) +2l~yr{y) \n[2{y + I)] + 2yh{y) - A + 0{t] 



;i8) 



The functions r{y) and hiy) are given by 



r{y) 

h{y) 



1 



L^il - yl) - L,{1 - y_)\ +-^Jy'-l r\y) 



(19) 



where y± = y ± — 1, and L2{x) is the dilogarithm (see Appendix B). They satisfy 
r(l) = h{l) = 1. 
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B. Gluons attached to the hght quark Hne 

Next consider the self-energy contribution for the hght quark shown in diagram D4 in 
Fig. |l]. It gives 

= -A{D - 2)N^CFgjTT{V'TV-f''-f'^-f"'} 

X fdSdi ^-^^^ 2 . (20) 



{00 + 2v s){uj' + 2v' ■ s){sYt^ {s - t) 



After performing the integral over t, the remaining one-loop integral can be readily evaluated 
using the master equation (|AT|). After Borel transformation, one finds the simple result 

D, = ^ = --. (21) 

D-4: 2e ^ ^ 



C. Gluons attached to both heavy and light quark lines 

Let us now turn to the most cumbersome part of the calculation, namely the loop cor- 
rection of the heavy-light vertices. The contribution of the diagram in Fig. |^ is 

D, = 8N,Cf gl Tr{ v'v V 7"^7^} 

J {uJ + 2v■s){uJ + 2v■t){u' + 2v'■s)sH^s-ty' ^ ' 

It is convenient to split the calculation into three parts by use of the trace identity 

Tr{P'rP7"^7'^} = Tr{p'rp((7"'3 + i[7",7'3] + 2t;V)}- (23) 

Let us denote the corresponding contributions by D^^ and discuss them in turn. 

1. Calculation of D^^ 

By rewriting its numerator, the integral appearing in dI^^ can be further decomposed 
into three parts: 

{u + 2vs){u + 2vt){u' + 2v'-s)sH^{s-tY = + + /s • (24) 

The integral Ji can be evaluated by introducing a new variable s' = s + t and using the 
master equation ([Ml). We find 



, / .\ 2 ^-'^ 

'^ = %mrni-o)m-i)l<iu-, 

^ ' 

X /dt ^- . r-r^. (25) 
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where V = {1 + 2yu + u'^Y^^, Q = lo + uu', and Va = {v + uv')a/V is a unit vector. The 
integral over t can again be performed using the master equation, resulting in a double 
parameter integral. The result simplifies upon Borel transformation. We find 

The calculation of the remaining integrals is straightforward. For I2 we introduce t' — t-\-s 
and make repeated use of the master equations. factorizes into the product of one-loop 
integrals. To evaluate its Borel transform it is convenient to combine denominators by a 
Feynman parameter. This gives 

f _ 2C 1 

I9 — 



{D-3){D-2) [2(y + l)]^/^-^' 

G(2-A0,f -l;y). (27) 



{D - 2) 

The calculation of the parameter integral in 73 is discussed in Appendix B. 

2. Calculation of D^^ 

For the two-loop integral appearing in D^^ we use the integral representations discussed 
in Appendix A. We find that 



iTr{7^'r7^[7",7^]} Jdsdi 



{u + 2v s){u + 2v ■ t){u' + 2v' ■ s) s'^t'^ {s - tf 
^{y-l)TT{V'TV}h, (28) 

where after Borel transformation 

C f^.o r.r. D ^ G(l,0,f;|/) 



G(3-A0,f;y)- 



D-2,\ ' ' '2'-^ [2{y + l)fl'-' 

r(D-i) 7^ 7^ - 1)^/'-' 1 



r(f)r(f-l)J V ^[u, + 2{y + l){u2-l] 

The double integral becomes trivial in the hmit D — > 4. The evaluation of the first parameter 
integral is outlined in Appendix B. 

Remarkable cancellations take place when one adds up the contributions from Ii to I4. 
We find the simple result 

+ = A{-i + 0(6)}. (30) 
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3. Calculation of D 



(3) 



This part of the amphtude involves the hardest integrah 

2v ■ t Sj3 



J'^''^\uj + 2vs){uj + 2vt){uj' + 2v' ■s)sH^{s-tY' ^^^^ 

It can be simphfied by using the method of integration by parts, which allows one to reduce 
a given loop integral to a series of simpler integrals |l9| , pO| , pTl . In this case, we can relate 
to a sum of four integrals involving five (instead of six) types of propagators. To this end, 
we evaluate the identity 



dsdt 



~ d 



2v ■ tsp {t 



dta {uj + 2v- s){uj + 2v ■ t){uj' + 2v' ■ s) s^t"^ {s - ty 







(32) 



to obtain 



(D-4)J^= dsdi 



S0 



- I dsdi 
dsdt 



{uj + 2v s) {uj' + 2v' ■ s) s2 12 (s _ t)2 
2v ■ t sp 

{uj + 2v s){uj + 2v ■ t){uj' + 2v' ■ s)(t2)2 (s _ t)2 
2v -tsp 

{uj + 2v s){uj + 2v ■ t){uj' + 2v' ■ s)(t2)2 s2 



{uj + 2v ty{uj' + 2v' ■ s) s2 12 (s _ t)2 



j{2) j{3) j(A) 

•^P '^P "T "^Z? "T ■-'13 



(33) 



The first three integrals can be calculated along the lines discussed above. After Borel 
transformation, we find 



J, 



(1) 



4C 



{D-A){D-2) [2{y + l)] 



D-2 ' 



J, 



(2) 



c 



[2(2/ + 1)] 



D/2-2 



{v + v')p 



2vr- 



D-2 



G(0,0,f -l;y) 



2(y + l) 

+ G'(0,l,f;y)t;;3 + G(0,0,f;i/)t;;,| 



J, 



(3) 



/3 



2C 



G(3 - D, 1, f ) t;^ + G(3 - D, 0, f ) t;' 



(34) 



2 7 ^/3J ' 

with C as defined in ([261) . Notice that because of the factor [D — 4) on the left-hand side of 
(^) these expressions have to be evaluated up to first order in e. 

The evaluation of j'jp is more involved. We have used the method of differential equations 
to calculate this integral ||22|. This interesting technique will be discussed below. For the 
moment we just present the result: 



J, 



(4) 



2Cv' 



13 



D-2\ I2(y + l)]^/2-i 
-2cjdu{l-u^-^^ ^"" + "')^ 



(D -4)G(0,D -4^ 



D 
2 



;i + 2yu + u^fl^ ■ 



(35) 
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We have not written the last integral in terms of G-functions in order to show explicitly that 
it is of order [D — 4). 

Next we set -D = 4 + 2e and expand the above expressions, keeping terms of order e. The 
integrals encountered are collected in Appendix B. They yield rather nontrivial functions of 
y. However, again remarkable cancellations appear if one sums the various contributions to 
the right-hand side of (PH]). Not only do the poles in 1/e contained in jj^^ and jj^^ cancel, 
but also most of the y-dependent terms. Our final result is, in fact, very simple. It reads 

C 



iD-A)if,TT{v'TV^^} 



[2(2/ +1)] 



X 4 - 2e 



(36) 



2 + 



27r2 



After Borel transformation, the diagram Dq gives the same contribution as D^. Hence 



A\--2 '-^ '^-^ 

e 3 



(y'-l)r\y) + 0{e)\. 



(37) 



D. Calculation of jj^^ using a differential equation 

In this paragraph we illustrate the application of differential equations to the analysis 
of multi-loop diagrams. Such techniques were introduced in Ref. ||2^ to evaluate integrals 
with massive propagators. They are readily adaptable to HQET. The idea is to derive 
a differential equation for a particular loop integral whose inhomogeneous term can be 
calculated in terms of simpler integrals. The original loop integral is then obtained from the 
solution of the differential equation. 

We start by rewriting jjp = x-^Ji3{x)\x=i, where 



Jsix) 



dsdt 



{xio + 2v ■ t){uj' + 2v' ■ s) s2 12 (s - ty 



(38) 



To derive a differential equation for J/s^x), we use again the method of integration by parts. 
Starting from the identity 



/ 



ds dt — 



dta (xu + 2vt) {iu' + 2v' ■ s) s2 12 {s-ty 







we find 



where 



ds dt 
dsdt 



Sf3 



(39) 



(40) 



{xuj + 2v ■ t){uj' + 2v' ■ s) s2 [(s - tyy 



{xuj + 2v ■ t){uj' + 2v' ■ s) t2 [(s - tff 



(41) 
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does indeed only contain simpler integrals. The general solution of the differential equation 
is 



Jr(x) = X 



<lzz^-''fp{z) + kp 



(42) 



with kjs being independent of x. We are interested in the Borel transform of this equation. 
A straightforward calculation gives 



(47r) 



D 



[1 + 2yz + z^f^^-' 
D-2 



(43) 



l-z 



A-D 



^D-3 _^ ^^/^^ 

;i + 2yz + ^2)^/2 



In order to determine the constant of integration we consider the limit x — > oo in ( p8|) 
in which 



lim xJji{x) = — Idsdt 



Sf3 



00 J {uj' + 2v' ■ s) {s - ty 



r(f-i)r(f-2).' 



\2D-5 



By evaluating (|42|) for x ^ 1, on the other hand, we find (for D > 3) 

lim xJp{x) = ^ ^/^^^ r(f - 1) r(f -2)v'^ + kfs Jim , 



(4vr; 



D 



Hence ki^ = follows. From the solution of the differential equation we now obtain 



(44) 



(45) 



d 



x=l 



-/^(l) + (D-4) I dz z'-^ f^iz) . 



(46) 



Substituting here z = 1/u leads to 



E. Summary and renormalization 

We are now in a position to sum up the various two-loop corrections that contribute at 
order as to the Borel-transformed correlator. We find 



6 

1=1 



A 



3/1 Att^ 



2\e 



9 



yr{y) - 1 



2 In [2(2/ + 1)] 



+2 



yh{y)-l -{y'-l)r\y) + 0{e) 



(47) 



Next we expand A from (|1J) around D = A, keeping terms of order e, and relate Si to the 
lowest-order correlator 
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This yields 



Si = — So -T+ 21n^ + 27. 



TT 



e 

7(y) 



47r2 \y + l 



Tr{ P r P } 



(48) 



T 



47r2 17 

~9~ ^ y 



-T + (21n| + 27^,-3 



+ Cpert(y) + C(e) 



(49) 



where 

1 

e 

We have introduced the functions 

4 
3 

i{y) 



1 , Att 

- + 7£; - In — 



(50) 



(51) 



^pcrt 



4 In 2 - 3 + In 



y + l 



4r 
3 



49 



-ln2 - 

9 54 



15 



In 2 



yh{y) - 1 
197 



+ In 



600 



iy-ir + .. 



both of which vanish at y = 1. The first two terms in the expansion of Cpertiv) around zero 
recoil were previously calculated in Ref. []TB[, and we confirm the result obtained there. 

The 1/e poles in (|49|) cancel upon renormalization of the heavy-light and heavy-heavy 
currents in (|l]). In the MS subtraction scheme, the corresponding renormalization factors 
are [|,|3|,i 



2i 



hi 



27re 



Z 



hh 



1 



as 



i{y) ■ 



(52) 



That means that our two-loop calculation reproduces correctly the known running of the 
hadronic form factors -F(/i) and ^(y, /i). By comparing (|49D with ([TT|) we can now write our 
final result for the renormalized correlator in form of a dispersion integral and introduce the 
continuum threshold uiq to obtain 



3 / 2 \ 2 7 



X 1 + 



TT 





/i 47r^ 17 / N , A* 

2 In + -^ + ^- IW In + Cpert(l/) 

9 3 6^;+ 



(53) 



This is the exact expression for the perturbative part of the correlator at order a^. It is now 
seen that at zero recoil one indeed recovers (p!0|). 
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IV. RENORMALIZATION-GROUP IMPROVEMENT 
AND NUMERICAL ANALYSIS 



The theoretical expression for the correlator depends on the subtraction scale /x, indicat- 
ing a scheme-dependence associated with the subtraction of the 1/e poles. This just reflects 
that the hadronic parameters F{fi) and ^(|/, fi), which are defined in terms of matrix elements 
of currents in the effective theory, are scheme-dependent quantities. At next-to-leading or- 
der in renormalization-group improved perturbation theory one can define renormalized, 
scheme-independent form factors by 



2/9 



1 - 



TT 



hi 



-a-Liy) 



1 - 



TT 



(54) 



where aL^y) = §7(1/) [H, and we have used that the number of light quark flavors in the 
effective theory is rij = 3. The next-to-leading logarithmic corrections consist of two parts. 
The coefficients Z are renormalization-group invariant quantities. For Uf = 3, they are given 
by @J2|J2|| 



185 7tt^ 



Zhhiy) 



324 243 ' 

7(y) [53 _ _ 
3 154 12 



(55) 



yr{y) + 



1 



L2il-y'_)+\n\y^ 



,,,r, if, ,2/1 smh 9 cosh 9 . smh 9 

dip ip coth yj — 1 itp coth 9 H ^ T~ — 1 — 

9 J L J I sinh 9 — sinh ijj sinh ip 



/752 Stt 
V729 ~ "8 



/ 368 47r'^\ 
V1215 ~ 135 



j(y-l)^ + ..., 



where y- = y — \/y'^ — 1, and the hyperbolic angle 9 is defined hj y = cosh 9. The coef- 
ficients 6, on the other hand, are scheme- dependent. They arise from matching of QCD 
onto the effective theory and combine with the scheme-dependent terms in (^) to give a 
renormalization-group invariant result. In the MS subtraction scheme one has 



4,26 



5i 



hi 



2 
3 



5hh(y)=0 



(56) 



After renormalization-group improvement, the sum rule for the renormalized Isgur-Wise 
function takes the form 



^ren(y) = «.(T) 



-aUy) K{T, uJo,y) 



(57) 



where we have summed the leading logarithms down to a characteristic scale given by the 
Borel parameter T. The renormalization-group invariant function K is given by 



13 



3T3 



TT^ \y + 1 



u)q/T 



X 1 + 



as{T) 



TT 



^ + y - 2Zhi - [2 - 7(2/)] Inx + Cpert(y) - ^hh(2/) 



a.(T) 



TT 



2Zhi + 7(?/)[Ei(-^)-7ij] +C(gg)(|/)-Zhh(|/) I 



^ ^ y-l \ {a,GG) ^ (2y + 1) (gg) 



.?/ + iy 487rT '3 4T2 ■ 

We have included the order-cts corrections to the quark condensate as calculated in Ref. ||12||.l 
The function C(^gg){y) has a similar form as Cpeniy)- It reads 



i{y) 



41n2 + ln 



y + l 



56 



V 9 27 



)(y-l)-( 



yh{y) 

-ln2 

15 225 



(y - ^)r{y) 



ln2-H|)(y-ir + 



(59) 



A remark is in order concerning the appearance of 7^; and the exponential integral in (pSD- 
The effective spectral density for the quark condensate receives contributions proportional 
to S{u!+) and l/u^ + S{uj^) ln(ci;+). The latter ones have to be regularized in the dispersion 
integral, leading to 



lim 



+ S{u!+ — rf) In 



In- 



T 



lE ■ 



(60) 



If one introduces a continuum threshold, one obtains an extra contribution Ei(— c^o/T) from 
the second term, where Ei(— x) = — ye~* is the exponential integral. This contribution 



is very small and has been neglected in Ref. |jT2[. 

For the numerical analysis of the sum rule (|57D we use the vacuum condensates as given 
in d^), as well as = 0.25 GeV (for = 3) in the running coupling agiT). In Fig. ^(a) we 
show the range of predictions for the renormalized Isgur-Wise function obtained by varying 
the continuum threshold over the range 2.0 < ujq < 2.6 GeV, and the Borel parameter inside 
the "sum rule window" 0.8 < T < 1.2 GeV, where the theoretical calculation is reliable. 
This window is determined by requiring that the nonperturbative contributions to the sum 
rule be less than 30% of the perturbative ones (T > 0.8 GeV), and that the pole contribution 
account for at least 30% of the perturbative part of the correlator (T < 1.2 GeV). The lower 
limit on T also ensures that as{T) is small enough to allow for a perturbative expansion. 
The above range of values for ujq was obtained from the study of the correlator of two 
heavy-light currents. As in previous analyses, we observe excellent stability of the sum rule. 
The sensitivity of the Isgur-Wise function to different choices of the continuum model is 
investigated in detail in Refs. |T2|,|T5[. We do not discuss this subject here, since we are 



The function ciq„\{y) was called ■nCy^{y) in Ref. |12|. Note that we have simplified the diloga- 



rithms appearing in h{y) as compared to this reference. 
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mainly interested in the effects of radiative corrections. We just note that the theoretical 
uncertainty in the sum rule prediction is probably larger than indicated by the width of the 
band in Fig. ||(a). 

To study the importance of radiative corrections we first have to relate ^rcn(z/) to a 
more "physical" form factor, which includes the logarithmic dependence on the heavy quark 
masses. Otherwise it is not possible to consider the limit — * 0. For simplicity, we work 
with a single scale m and define 

Cphys(?^, y) = [a,(m)] ''^^^jl + ^^il^l Zhh(l/)| ^ren(l/) • (61) 

We use m ~ 2.3 GeV as a characteristic scale for 6 c transitions Fig. ^(b) shows the 



next-to-leading order result for this form factor in comparison with the "bare" Isgur-Wise 
function computed by neglecting radiative corrections. We also show the leading logarith- 
mic approximation to ^phys{^,y), which is obtained by ignoring terms of order in (|58D 
and keeping only the first factor in (|6TD. It is apparent from this figure that the radiative 
corrections to the Isgur-Wise function are well under control. The large y-independent cor- 
rections, which enhance the sum rule prediction for the decay constant F by 50% [jTT],|T2|,|13[ , 
cancel out in the ratio (|57|) . The remaining recoil-dependent radiative corrections are small. 
At large recoil, they tend to decrease the form factor by 5 — 10%. Part of this effect comes 
from leading logarithms and is associated with the velocity-dependent anomalous dimension 
aiiy) of the heavy- heavy current in the effective theory. 



V. CONCLUSIONS 

We have presented the complete QCD sum rule analysis of the Isgur-Wise form factor 
^(f ■ v') at next-to- leading order in renormalizat ion-group improved perturbation theory. 
To this end, we have derived the exact result for the two-loop corrections to the triangle 
quark loop. Such a calculation, which was never done before for a form factor of heavy 
mesons, becomes feasible by using the heavy quark effective theory. We have developed 
some general techniques for dealing with two-loop integrals with two different types of heavy 
quark propagators. Using the method of integration by parts, complicated integrals can be 
reduced to simpler ones in a recursive way. Integrals which cannot be reduced any further 
can be evaluated by using differential equations. We have applied this technique to the loop 
corrections to the heavy-light vertices. We have also presented an integral representation for 
two-loop integrals which is particularly convenient for QCD sum rule calculations. These 
methods can be applied to other sum rule calculations and will eventually lead to more 
accurate predictions for heavy meson form factors than were available before. 

Our numerical analysis shows that, unlike in the case of meson decay constants, radiative 
corrections to the Isgur-Wise function are small and well under control. This is an important 
result which puts the sum rule analysis of ^{v ■ v') on a firm footing. The smallness of 
the two-loop corrections in this particular case was not unexpected, however, since the 
normalization of the Isgur-Wise function at zero recoil prohibits any recoil-independent 
radiative effects. This does not imply that such corrections are always negligible. In fact, 
some of the universal functions appearing at order l/mg in the heavy quark expansion 
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receive their leading contributions at order a^. Then the two-loop perturbative contribution 
is important and cannot be neglected for a reliable analysis . 
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APPENDIX A: LOOP INTEGRALS IN HQET 



1. One-loop integrals 

We summarize some important equations for one-loop tensor integrals in HQET. Integrals 
involving massless propagators only can be found, i.e., in Ref. p8[. Integrals involving two 
types of heavy quark propagators were considered in the second reference in . There the 
master equation 



//.l...Mn(",/5,7) 



tfll ■ ■ ■ tfj.„ 



(Al) 



/„(«,/?, 7) / du 



u 



7-1 



1/3+7 



-I D-2a+n 



V{u) 



(47r)^/2 

was derived, where dt = {27!-)^^d^t, and 

r(2a + P + -f-D-n) T{D/2 - a + n) 



^W...Mn(M) , 



/„(«,/?, 7) 



r(a)r(/5)r(7) 

Q{u) = UJ + uuj' , 

V{u) = {l + u"^ + 2uvv'y^^ . 

For n = 0, 1, 2 the tensors K^^,,,^^{u) are given by 

K{u) = 1 , 



K,,,{u) = V^{u)V,{u) 



9 



D~2a + 2 



with Vi_i{u) = {v + uv')^/V{u) being a unit vector. We note that the master equation is valid 
for arbitrary values of a, /3, and 7. 

In the case of one heavy quark the master equation reduces to M, 



•■fll-.-l-lr. 



dt tfj,^ . . . tfj,^ 



UJ 



-f^J \uj + 2vt 



D-2a+n 



where 



(47r)«/2 

_ r(2a + p-D-n) T{D/2 -a + n 

^ r(a)r(/?) 



(A2) 



K^^,,,^^ is obtained from above by replacing V^{u) by v^. 

The integral representation (XT) is particularly convenient for a Borel transformation in 
UJ and uj' . Defining the Borel operator by 



17 



we note that 



2. Two-loop integrals 

Let us now turn to two-loop integrals in HQET. The case with one heavy quark has been 
discussed in detail in Ref. ||2^. Here we consider integrals with two species of heavy quark 



propagators. They have the general form 

C::Z("'/^'7'^;«'^'C)= (A3) 

^fll ■ ■ ■ O/lm ^ • . . t 



ds dt' 



{uj + 2v ■ sT{uj + 2v ■ tY{uj'+2v' ■ sy{uj' + 2v' ■ ty {-s^)\~t^y[-{s - tyf 

We will derive a representation for this integral which is particularly convenient for further 
analysis. The first step consists in performing a Wick rotation of the loop momenta, 



s 



t ^ ^t\t), 



so that 



2 / 0\2 I ->2 2 

V ■ s iv^s^ — V ■ s = we-se , 

where se = (s", s ) and ve = iy^.iv) are vectors in a Euclidean space. Note that this 
definition of a Euclidean velocity ensures that v% = 1 and ve ■ v'^ = y, where y = v ■ v' in 
Minkowski space. After the Wick rotation, we represent the massless propagators as Fourier 
integrals in a Z)- dimensional Euclidean space and use an exponential integral representation 
for the heavy quark propagators: 



1 i-iy 



{uj + 2ivE-SEy ri 



a] 







The most general two-loop integral involves three £)-dimensional integrations over Xi and 
four one-dimensional integrations over Aj. The advantage of these representations is that 
the integrals over the loop momenta can immediately be performed and give rise to two 
D-dimensional 5-functions, which eliminate two of the integrations over Xi. Furthermore, 
note that two of the integrals over Aj become trivial upon Borel transformation, since 
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The tensor structure in the numerator in can be generated by taking derivatives with 
respect to Xj. Recall that for every timelike index there is a factor i encountered during the 
Wick rotation. For every spacelike index, on the other hand, one encounters a factor i when 
rotating back to Minkowski space. Together with a factor from the loop integrations there 
is thus a factor of taken into account. 

Let us illustrate this technique for some of the integrals encountered in Sec. |T|. We 
start with the integral in (0). It contains only two types of massless propagators, and 
consequently there are no integrations over Xi left after evaluation the (5-functions arising 
from the loop integrations. In the notation of ( [A3| ) we find 

Ul, 1, 1, 1; 1, 0, 1) = ^^^\2t~ Jd^i dAs dA4 e(^^+^^)-+(^3+A4K 



(4vr)^ 

Xia 1 

X 



where Xi = (Ai + \2)v + (A3 + A4)f and X2 = X2V + X^v'. After Borel transformation we 
obtain 

4(1, 1, 1, 1; 1, 0, 1) = ^^^''^''1° / d^i dz2 ^ 



with C as defined in (EB). Substituting Z2 = uzi, an integration by parts in zi yields 



2 



/I z 
d^i -jjj^ = . G(0, 0, f - 1; y) , 

Q [zi + zi + 2yziZ2) ' l^-4j 

which leads to (|TB|). The integral Jg in (^) can be evaluated along the same lines. In this 
case there are only three integrations over Aj, and one is thus left with a single parameter 
integral after Borel transformation. 

A more complicated integral is that appearing in (^). Following the general procedure 
outlined above we derive 

/|(1, 1, 1, 0; 1, 1, 1) = ^ ^ / dAi dA2 dAs e(^^+^^)-+^3^' 



(47r)^/' 







x\fl' {xlf/' [xlf"-' ' 

where X2 = Xi + (Ai + X2)v + X^v' , and X3 = Xi + Aif + Asf'. The integral over xi has 
the form of a Euclidean one-loop integral and can be performed in the standard manner by 
introduction of two Feynman parameters Zi. One then contracts the Lorentz indices with 
those in the trace in ( p8D to compute the integral 1^. The result is 



J dX, dX2 dAs A.As e J dz, dz2 -J^^^^^^J^ > 
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where Zi = 1 — Zi, and 

M'^{Xi, Zi) = Zi{z2p'^ + Z2q^) - zf{z2P + Z2qf 

with p = X2V and q = (Ai + X2)v + X3V'. After Borel transformation the integral can be cast 
into the following form: 



T{D - 1) JdX X^-^ J dui J du2 



For D < 3 one can use an integration by parts in A to obtain (ppj). By analytic continuation, 
this result can then be evaluated around D = 4. 



APPENDIX B: PARAMETER INTEGRALS 

We collect some useful formulae for the evaluation of parameter integrals. We start with 
a remark on divergent integrals such as those appearing in (^) and (^). Assuming first 
that D is sufficiently small, one can use an integration by parts to rewrite these in terms of 
integrals which have a well defined expansion around 0=4. For instance, for D < 4 one 
can show that 



{D-4)G{3-D,0,^;y) = D\yG{4-D,0,^ + l;y) + G{4-D,l,^ + l;y) 



1. 



Similarly, G(2 — 0, y — 1; y) can be related to G(3 — D, 0, y; y) plus nonsingular terms 
for D < 3. By analytic continuation, one can then evaluate the resulting expressions in the 
vicinity of D = 4. 

We now present a list of parameter integrals which are encountered when one expands 
the results presented in Sec. fT| around D = 4. When evaluating these integrals it is useful 
to introduce a hyperbolic angle 6 hj y = cosh 6. Then 

V^{u) = l + 2yu + u'^ = {u + e^){u + e"^) 

factorizes. Setting Rq = = 2{y + 1), we find: 

1 

1 r{y) 



Fi= Jdu 



, V\u) 2 ' 
Ro jdu ^ , = \"^^'! +l. 





1 





1 



Ro du = _!M _ Ur(y) + 1] lni?o -2l{y)} 
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du 



r{y) 



\nRo + h{y) 



Rn / du 



[V^{u)f ~ 2{y - 1) 



F7 = du 



ln(l -n) 



Rl 



{2y - 2h{y) + [r{y) - 1] (1 + lni?o)} , 



In^i^o y^-l^,,_, y + 1 



ylnRo {2-y){y+l) 



2(y-l) 4(y-l) 
\nV\u)-\nRo 



8 12 
r{y) \nRo-2l{y) 



■r (y) 



r{y) 



R^Q Idu 



[l-u) [V^{u)r 



6 



4 



r^{y) 



y + 1 



r{y) \nRo 



+ {y + 1) (1 - lni?o) - - + + 1) [F5 + Fe - lni?o 



The functions r{y) and h{ii) have been defined in (p^. In addition, we have introduced 

1 



Ky) 



which satisfies Z(l) = 2 In 2. Here L2{x) = — /q yln(l — t) is the dilogarithm. The first 
derivatives of these functions at y = 1 are r'(l) = — |, /i'(l) = ^, and /'(I) = | — |ln2. 
Finally, we note the following useful identity {n > 1): 



1 

y du ln"(l 



_D(y + M) 



D/2+1 



n du 



In^-^l -m) 



-D/2 
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FIGURES 



FIG. 1. Two- loop diagrams contributing at order to the perturbative part of the sum rule 
for the Isgur-Wise form factor. Heavy quark propagators are drawn as double lines, while the wavy 
line represents the weak current. 

FIG. 2. (a) Sum rule prediction for the renormalized Isgur-Wise function ^ren(y)- The width 
of the band arises from variation of uq and T as specified in the text, (b) The "physical" form fac- 
tor '^phys(^)2/) computed in next-to-leading order in rcnormalization-group improved perturbation 
theory (solid), in leading logarithmic approximation (dashed), and without including any QCD 
corrections (dotted). We use the central values uoq = 2.3 GeV and T = 1.0 GeV. 
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